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Yukawa theory at vanishing temperature provides (one of the ingredients for) an effective descrip- 
tion of the thermodynamics of a variety of cold and dense fermionic systems. We study the role 
of masses and the renormalization group flow in the calculation of the equation of state up to two 
loops within the MS scheme. Two-loop integrals are computed analytically for arbitrary fermion 
and scalar masses, and expressed in terms of well-known special functions. The dependence of the 
renormalization group flow on the number of fermion flavors is also discussed. 

I. INTRODUCTION 

Cold and dense fermionic systems are found in a wide range of physical environments and energy scales. In condensed 
matter, for instance, in the phenomena of antiferromagnetic ordering and superconductivity in the Hubbard model 
Increasing dramatically the density, Fermi pressure is responsible for compensating gravity in the hydrostatic 
equilibrium of compact stars, the Fermi gas being formed by electrons, neutrons or even quarks, depending on 
the density and formation process Q. In nuclear and particle physics, hadrons undergo chiral and deconfinement 
transitions at sufficiently large baryonic densities 0], exhibiting a very rich phenomenology including several color 
superconducting phases according to the different possibilities for quark pairing [i| . 

Unfortunately, the thermodynamics of such plasmas in the case of vanishing temperature and finite density is not 
very amenable to first-principle calculations. On one hand, the so-called sign problem brings about major technical 
difficulties for performing Monte Carlo lattice simulations at finite chemical potential [f|. On the other hand, although 
the perturbative series for the thermodynamic potential at zero temperature and finite density seems to be much 
better behaved than its counterpart at finite temperature 0, @, HQ, the values of chemical potential that are 
phenomenologically interesting usually have some overlap with a region where the coupling becomes large due to 
renormalization group running, so that perturbative calculations break down. Therefore, one is usually compelled 
to resort to effective theories, either to simplify the original description in terms of a fundamental theory, but still 
keeping its relevant symmetries, or to complement a given perturbative approach in the region where the coupling 
becomes large. 

In spite of the fact that an effective theory does not require renormalizability to be well posed, this attribute 
is highly desirable if one is interested in investigating the behavior of physical quantities as the energy scale is 
modified, as discussed above, which can be accomplished via renormalization group methods. The requirement of 
renormalizability restricts considerably the spectrum of possibilities, and in most cases the effective theory will contain 
a Yukawa sector. In fact, the Yukawa theory is very common in particle and nuclear physics, not only because it 
is very convenient in building effective models but also due to the need for a mechanism of spontaneous symmetry 
breaking and mass generation in gauge theories. In condensed matter systems exhibiting phase transitions, such 
as the antifcrromagnetic-supcrconductor transition mentioned above, the situation is analogous. Thus, perturbative 
Yukawa theory in a medium has a wide range of applications. When it represents a sector of a fundamental theory, 
such as the quark-Higgs sector of the electroweak theory, one must assure that the scales are such that the Landau 
pole is not reached within the region of investigation, so that one can compute loop corrections and implement a 
renormalization group (RG) analysis. On the other hand, when it plays the role of part of an effective theory, one 
can avoid the Landau pole by simply introducing a cutoff beyond which the theory is meaningless. In the latter, one 
can choose to take effects from RG flow into account or not, depending on the system under investigation. 

In either case, the perturbative treatment of the Yukawa theory at finite density may be of help also for testing 
lattice simulations in the limit of small coupling. In fact, numerical studies at zero temperature and density of the 
Yukawa theory with a real scalar field [l(| and of the chirally invariant Higgs- Yukawa model have revealed a rich 
structure, in spite of the formal simplicity of the theory. In some cases perturbative calculations may complement in 
an efficient way Monte Carlo simulations in the investigation of the phase structure of a given effective theory fl2| , in 
others the effective model can use lattice data as input for its parameters 
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In this paper wc investigate the role of masses and the RG flow in the calculation of the equation of state of cold and 
dense Yukawa theory up to two loops within the MS scheme. Previously we have presented preliminary results on the 
two-loop correction for the pressure with one massive fermion flavor and a massless scalar field [3] and the influence 
of the RG running on the equation of state for Np flavors [l5| . The interest in studying results for different numbers 
of flavors, such as Np = 4, is motivated not only by the phenomenological interest in different physical systems 
but also by on-going studies on the lattice using Kogut-Susskind fermions [l6| , Here we present full results for the 
thermodynamic potential with arbitrary masses and number of fermion species, including the effects from RG flow of 
couplings and masses. The two-loop momentum integrals are computed analytically for arbitrary fermion and scalar 
masses, the final result being expressed in terms of well-known special functions. Renormalization is implemented in 
the standard fashion in the MS scheme. 

The paper is organized as follows. In Section II we show our results for the two- loop expansion of the thermodynamic 
potential of the Yukawa theory at finite density, and discuss the renormalization procedure. Fixed mass and coupling 
results are presented in Section III. In Section IV we study the coupling and mass RG flows and their influence on 
the pressure. Section V contains our conclusions and outlook. Several technical details are left for a final appendix. 



II. TWO-LOOP THERMODYNAMIC POTENTIAL 



In what follows, we consider a gas of Np flavors of massive spin- 1/2 fermions whose interaction is mediated by a 
massive real scalar field, 0, with an interaction term of the Yukawa type, so that the Lagrangian has the following 
general form: 



Cy — + + C int , (1) 



where 



= ^ot - m ) V'a , (2) 

£<f> = i(^)(9^)-im^ 2 -A 3 3 -A0 4 , (3) 

N F 

C mt = tp a tp a (j). (4) 

a=l 

Here, m and m<f, are the fermion and boson masses, respectively, assuming all the fermions have the same mass, for 
simplicity. The Yukawa coupling is represented by g; A3 and A are bosonic self-couplings allowed by renormalizability. 
The latter play no role in the thermodynamics up to this order, unless in the presence of a nonzero scalar condensate 
[l7| . as in the case of a spontaneously broken symmetry, where the condensate contributes to the effective masses. 
In this paper we assume (cf>) =0, leaving the treatment of spontaneous symmetry breaking, as in the linear sigma 
model, for a future publication [18j . 

Although we are interested in the limit in which the temperature T vanishes, it is technically more convenient to 
work with the imaginary time formalism of finite-temperature field theory, where the time dimension is compactificd 
and associated with the inverse temperature (3 = l/T [ljj]. At the end, we shall take the limit T — * 0. As is well-known, 
to be consistent with the spin-statistics theorem for bosons (B) and fermions (F), one has to impose, respectively, the 
fields to be periodic or anti-periodic in the imaginary time r, so that only specific discrete Fourier modes are allowed. 
Therefore, as is customary in finite-temperature field theory, integrals over the zeroth four-momentum component 
are replaced by discrete sums over Matsubara frequencies, denoted by ui^ = 2nirT and u>^ = (2n + 1)ttT, with 
n integer. Taking finite density effects into account amounts to incorporating the constraint of conservation of the 
fermion number which, in practice, is implemented by a shift in the zeroth component of the fermionic four-momentum 
p° = iujn *-* P° = i^n + I 1 ' P" being the chemical potential. 

From the partition function written in terms of the euclidean action for Eq. ((T|), Zy(T,fi) = Tr cxp(— Sy), one 
derives the pcrturbative series for the thermodynamic potential Qy = —(l/j3V) \nZy: 



Civ = In Zn In 

PV (3V 



(5) 



where V is the volume of the system, Z is the partition function of the free theory and Si nt represents the euclidean 
interaction action. Notice that Wick's theorem implies that only even powers in the above expansion survive, yielding 
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a power series in ay = ff 2 /47r [20I ]. Omitting purely bosonic contributions (which are /^-independent) and the diagrams 
representing counterterms, the thermodynamic potential up to two loops is given, diagrammatically, by 



1 

W 



N F 



o 



(6) 



where solid lines represent the fermions and dashed lines stand for the bosons. The first diagram corresponds to the 
free gas contribution 



o 



-2V 



(2tt) 3 



PE p + In (l + e -^p-^)) + hi (l + e ~^ E ^ 



(7) 



with 



1/2 

<p 2 + m 2 Y'. The zero-point energy divergent term can be absorbed by a convenient redefinition of 



the zero of the thermodynamic potential, since it is independent of T and /i. The O(ay) correction is given by the 
exchange term (l9j : 



e = (3V g 2 4 Tr 
Tp u p 2% k 



{2ir) 3 f3S^(K-P 1 +P 2 ) 
f> 1 -m){ml-K*){f 2 -m) 



(8) 



where the trace is performed over the Dirac structure, and the 4-momenta are given in terms of the Matsubara 
frequencies and the 3-momenta for fermions, Pi = (p° = iu^. + /i , pj) , and bosons, K = (fc° = iuf , k) . We use the 
metric tensor g^ 11 = diag(+, — , — , — ) and the following notation for the sum-integrals: 



(9) 



The calculation of the exchange diagram is very similar for different theories and is a standard exercise in finite- 
temperature field theory, the main difference in details coming from the tensor structure of each theory under inves- 
tigation. Previous results for cold and dense systems can be found, for instance, in Ref. 1211 for nuclear matter, in 
Refs. 0, [13, [H for QED and in Refs. 0, M, HE E3, H HI for QCD (see also Ref. [Hj). Aft er some long but 
straightforward algebra, we can write the exchange term for the Yukawa theory in the following form (see Appendix) : 



e 



1 



rf 3 Pid 3 p 2 

(2tt)6 2uj 12 E p1 E p 



J + ^12 Si + J - W12 £2 



J- E+ - J+ E_ 



n b (cu 12 ) N f {l) 



- J + (E-+lo 12 )- J-(E + -uj 12 )\ N f (l) 
-2 J_ E+ n b (cu 12 ) - 

-J_(£+-w 12 )j , 

where we have defined the following functions: 

o m 2 - Pl -p 2 ±E pl E P2 4m 2 - 



J± 

N f (i) 
Si 
£2 



, ,2 

E T - uj 12 



E^ - u> u 



n f (E Pl + (j,) n f (E P2 + /i) + n f (E pl - //) n f (E P2 
n s {E pl + fi) n s {E P2 - fi) + n f (E pl - fx) n f (E P2 

1/2 



m) 
m) 



with E± = E P1 ± E P2 and uj 12 = (Jpi - p 2 | 2 + m 2 A ; n b {uj) = [exp(/?cj) - 1] and n s (E) = [1 + exp(f3E)}~ 



(10) 



(11) 

(12) 

(13) 
(14) 

1 

are 



the Bose-Einstein and the Fermi-Dirac distributions, respectively. 
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The physical meaning of each term in Eq. ([10]) is clear. The first two lines arc quadratic in the statistic distributions, 
representing contributions coming from the scattering of particles from the medium. These are therefore ultraviolet 
finite due to the exponential suppression of the integrands implemented by the distributions. The other terms 
contain contributions coming from the scattering of virtual particles. The last term is a pure vacuum contribution, 
independent of T and /i, and can be absorbed by a redefinition of the zero of the thermodynamic potential. On the 
other hand, the remaining terms mix medium and vacuum particles, being linear in the statistic distributions, and 
call for rcnormalization. 

The divergent contributions that arc linear in rib, L^, and linear in rif, Ly, can be written in terms of amputated 
(amp) vacuum (vac) self-energies evaluated on the mass shell (m.s.) as follows (see Appendix): 

K 



-2 pV 



->b = 




(15) 



(16) 



Here, matter means that the pure vacuum part has already been subtracted, and, while K = (u>k, k) and 
Pi = (£ Pl ,Pi) are evaluated on the mass shell, Q = (iujf , k) and P = (iu;^ + /i,pi) are not. Implementing the 
rcnormalization procedure for the self-energies above in the MS scheme, one obtains the renormalized expressions (see 
Appendix) 



t ren 



t ren 



-2 PV^f- 4m 2 

47T 



cf' 



Pi 



(2tt) 3 



Nf(l) 



2E„ 



ml 



-2 pv 



ay 

47T 



d 3 k 

(2tt) 3 



2n b (w k ) 



2w k 



[2 a 2 + 6 a 3 ] 



(17) 
(18) 



where 



4m 2 
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+ tan -1 
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4m 2 



2m 2 



3 , (m 2 
2 V A^" 




o 2 



Q 3 



2^2 

~Q m 4 



2m 2 - 



12 



4m 2 



- 1 



4m 2 



tan 




4m 2 



4m 2 



4 m 2 



1 



log 



m 
A 2 



(19) 
(20) 

(21) 



and A is the renormalization subtraction point. 

In the limit of vanishing temperature and in the absence of a scalar condensate, the Bosc-Einstcin distribution 
vanishes, so that the purely bosonic diagrams, omitted in Eq. ((SJ), do not contribute. The Fermi-Dirac distribution 
simplifies to the Heaviside step function, limT->o n f(Pp ~ I- 1 ) = — ^p); signaling the occupation of all states in the 
interior of the Fermi surface. Simplifying all the terms in the thermodynamic potential, we obtain in the cold and 
dense limit 



-N F 



1 



24tt 2 



[2 MP/ 



3m 2 ul 



1 

— Nf ^TTCtY 



Jl 



1 



16tt 4 



m uct\ 



(22) 



where u = fipf — m 2 log j , Pf is the Fermi momentum pf = \/ [i 2 — m 2 and we have defined the integral 

d 3 p 1 d 3 p 2 1 



Ji 



(2nf 2E P1 E P2 



J+ 9(ji-E pl ) 9(ji-E P2 ). 



(23) 
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The evaluation of the integral above when both fields arc massive is highly non-trivial, and results are usually presented 
in the limit of at least one vanishing mass, which is the proper case in gauge theories, or numerically [l9| . We obtained 
the following complete analytic result for arbitrary values of m, and fi (see Appendix for details) 



2m 2 1 



4m 2 



T 2 

It U 

2 



(24) 



where the function Xj is given by 



2m 2 



log 



Ap 2 + m\ 



2 „2 



m % \ M Pf 



2m 2 



m 4> 
m 



4ra 2 



(fj, Pf + U) Aan +Pf — 2m 2 Klog 



Pf 



fi + m ' 4m 2 



(25) 



with 



Aan = tan 



2m(ji + m)y 1 



tan 



Pf 





4 m- 



. m<* \/ 1 



2m(/i + m) 



For x = pf /(fi + m) < 1 and z = m 2 i /4m 2 , we can write the function /Ci og (x, z) as 



/Ci og (x,z) = A(x, z) + -a/ z(z - 1) C L j(x, z) , 



(26) 



(27) 



where 

C u (x,z) 



-Li 2 



1-x 

1 + X 



Li 2 (l - [z - y/z(z-l)j (1 - x)J - Li 2 (l - [z + y/z{z-l)j (1 - x)J + 

l-x r 



1 - (z + Vz(z-l)) (1 - x) 



Li 2 



1 + X 



1 - (z - y/z(z- 1)) (1 - x) 









+ 


X I— i 


—x 



,(28) 



with the Dilogarithm Function [30( defined as Li 2 (z) = 'Yl^ =1 {z n / n 2 ) . The quantity A is defined in two regions, 
according to the ratio z, as 



A(x,z) 



A< (x, z) , if z < 1 
A>(x, z) , if z > 1 ' 



(29) 



where 



A<(x,z) = Vl - z{ Iog(l-x) 



tan 



-i(VT 



tan 



+ log(l + x) 



tan 



— tan 



Vz(l-z) 



A >(x,z) = -Vi - T< log 



' 1+x 2 
. 22 

r,2 1 



1 



1 — 6 [ mz 



(1 + x) 2 2mx 



l-x 2 l-x" 



(30) 



\ bg ( \y/z~(l + x 2 )-2xVz~T\ 



(1 + x) 2 J ° \ \^{\ + x 2 )^2x^z~T\ 

X 



[z(l - x 2 ) - 1] - 6 1 + 



(l-x 2 )7J(I^T) 



[2-z(l-x 2 )] + 



(l-x 2 )v^(^T) 



[2-2(1 



(31) 



With this general result for the thermodynamic potential for the cold and dense Yukawa theory, we can analyze the 
effects from the different nonzero masses and their competition, the role of the interaction and the renormalization 
scale, as well as the influence of the RG flow of the coupling and masses. 
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III. RESULTS AT A FIXED ENERGY SCALE 



Let us start our study by the simpler case with fixed energy scale A. In practice, one can choose an appropriate 
value for A, in this case a parameter in a given effective theory, by imposing, e.g., known experimental constraints for 
the system under consideration. In the next section, we investigate the role played by A as the running scale in the 
RG flow of g, m and m<$,. For fixed A, the main issue is the influence of the masses which we illustrate separating the 
analysis into two cases, the first with = 0, and then the general one. 



A. Massless boson 



When ?Ti0 = 0, the integrals involved in the computation of the thermodynamic potential at finite density simplify 
dramatically. In fact, the former complicated functions a\ and J\ reduce to 

„% Q ' = M log (S)- 1321 



m^— >o 32ir 

so that the thermodynamic potential assumes the much simpler form 



Sly = -N F ^[2w}-3m 2 u]-N F ^l { 3u 2 -4p} + m 2 u [ 7 -31og(^)]} 



(34) 



It is interesting to notice that the inclusion of mass for the fermions brings the presence of logarithmic corrections one 
order down in ay. This is a general feature, also manifest in theories such as QCD. In massless QCD, for instance, 
one has ~ a s corrections at two loops and ~ {a 2 s , a 2 loga s , a 2 log(A//x)} at three loops [1,0]. However, in the massive 
case one finds not only ~ a s terms, but also a contribution ~ m 2 a s log(A/m) at two loops [2{|, analogous to the one 
that can be seen above. In our case, this feature will be important when we incorporate the RG running of ay and m. 
On the other hand, since the explicit dependence on A is logarithmic, the effect of its variation without the RG flow 
does not affect significantly the pressure [l4| unless for very large values of the coupling, where perturbation theory 
is meaningless. 

Figs. [T] and [21 illustrate, respectively, interaction and fermion mass effects on the two-loop pressure, P = —Sly, as 
a function of the chemical potential \i in units of the reference scale A = A^^g. Fig. [1] shows that the pressure is 

raised as the system considered interacts more strongly, reaching a variation of ~ 50% in comparison with the free 
case for ay ~ 0.72 and (/x — m) = 0.1 Aj^g. This behavior observed at small chemical potentials is not maintained 
at large scales: there is a crossing at (fi — m) « 0.25 A^jg so that the interaction effects tend to reduce the pressure 

of highly dense media, although the corrections in this regime appear to be less significant 20%). In fact, since 
the Yukawa coupling is treated perturbativcly, it is quite reasonable that interaction effects do not modify drastically 
the thermodynamics within the domain of validity of our calculations. Nevertheless, depending on the specificities of 
the system under investigation, corrections of a few percent as the ones obtained above may have impact on physical 
predictions. 

On the other hand, the fermion mass is not linked directly to the expansion parameter and its effects are therefore 
not expected to be constrained by the perturbative framework. Indeed, Fig. [5] shows that the influence of finite 
fermion masses on the thermodynamics can be sensibly more consequential: at fi = 0.6A^jg, for instance, the 
inclusion of fermions with mass m = 0.4A-^g lowers the pressure by a factor of ~ 1/5 in comparison with the massless 
case. Even for masses one order of magnitude smaller than the reference scale A-^jg, the corrections are sizable for 
sufficiently small chemical potentials. Therefore, finite fermion masses alter significantly the thermodynamics of the 
Yukawa theory in a wide range of the parameter space and approximations that neglect them should be implemented 
cautiously. 

This result represents another indication of the potential importance of fermion mass effects. Recently, the modi- 
fications brought about by finite masses have received increasing attention in different contexts, either because they 
have been underestimated before or due to the interest in high-precision tests, in experiments and in more realistic 
lattice simulations. In calculations within the Standard Model, for instance, there is a whole literature on multi-loop 
integrals that is progressively turning its attention to the inclusion of contributions due to non-vanishing fermion 
masses |3lj . In QCD it was also pointed out recently that finite quark masses, especially the strange quark mass that 




FIG. 1: Pressure normalized by the scale A = A-j^g- 
as a function of the fermion chemical potential for 
different values of the coupling g. The fermion mass 
is fixed at m — O.IAt-th-. 



FIG. 2: Pressure normalized by the scale A = Aj^jt 
as a function of the fermion chemical potential for 
different values of the fermion mass m. The coupling 
is fixed at ay = 1/4-7T. 



is not that small compared to typical scales in QCD, should play an important role in the critical region of the chiral 
and the deconfining phase transitions Q , possibly bringing relevant astrophysical consequences [29| ■ These results for 
the Yukawa theory signal that also within effective theories mass effects bring relevant corrections. 



B. General Case 



Let us now consider the full massive case. The influence of the boson mass on the thermodynamics of the 
Yukawa theory is shown in Fig. [3J Since we are analyzing the cold and dense regime, the pressure investigated 
consists essentially in a quantum Fermi pressure. In this vein, it is convenient to interpret the results in terms of a 
free quasi-particle theory: the variation of the pressure can be seen as a consequence of the modification of the effective 
mass of the quanta present in the system through the radiative corrections to the fermionic self-energy due to the 
coupling with the bosons. We showed above, in Fig. [2 that for a given chemical potential the pressure for a theory 
of heavy fermions is lower than the one for light fermions. Therefore, Fig. [3] indicates that for m^/m sufficiently 
small the self-energy is negative, diminishing the fermionic effective mass and increasing the pressure. For the case 
illustrated in Fig. [3J this effect is inverted for m^/m w 15. As one raises the chemical potential this behavior is 
intensified, as shown in Fig. [4] 




However, it is important to stress that, from a quantitative point of view, these results depend strongly on the 
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rcnormalization scale A adopted. Although when analyzing a general theory it is not possible to discuss quantitatively 
the result, in the description of a specific system the renormalization scale is determined through physical constraints, 
such as positivity of the energy density and the domain of energies being investigated. In the case of QCD [f|, for 
example, it is not reasonable for quark matter to have a higher pressure than a hadron gas for low chemical potentials 
(i.e., equivalent to densities of the order of the nuclear saturation density, no ~ 0.16 fnG 3 ), since this would imply the 
existence of stable deconfined matter in a density regime where only hadrons and nuclei are observed. Equivalently, 
one expects the quark pressure to be higher than the hadronic one for densities w > In this way, one fixes a 
physically reasonable range for the renormalization scale in QCD. This kind of procedure can also be implemented 
once the phenomenology described by an effective theory containing a Yukawa sector is known, and the parameters 
are given by experimental observations. Therefore, once the physical picture is fixed, the predictions above can be 
precised quantitatively. 



IV. RG RUNNING EFFECTS 



The thermodynamic potential £ly depends on the rcnormalization scale A not only explicitly but also implicitly, 
through the RG running of the coupling and the masses, with A corresponding to typical momenta involved in 
scattering processes in the medium [32(. In the context of effective field theories, it is plausible to use the results 
above, with fixed ay and masses, and A determined by phenomcnological constraints. Nevertheless, for several 
applications, RG flow may be considered a relevant feature or even be intrinsically present, as in the case of sectors of 
the Standard Model. Then, one has to solve the corresponding RG equations and include their effects in the evaluation 
of the equation of state. In fact, these effects bring major consequences in cold and dense QCD [f| [2{|. Since we 
consider the general case of massive fermions and massive bosons, there will be an intricate competition between the 
different effects. 




FIG. 5: Yukawa coupling RG flow normalized by the 
reference scale Aj^g f° r different numbers of fermion 
flavors. 



FIG. 6: Fermion mass RG flow normalized by the 
reference scale Y^<t for different numbers of fermion 
flavors. We chose mo = O.lAj^g-, and vertical lines 
represent the scale at which ay = 1 for each value of 
Nf, increasing from left to right. 



In what follows we include the effects from running masses and coupling, and also investigate the role of the number 
of fermion flavors Np. Although we study the behavior of RG running for different numbers of flavors, our results 
for the pressure are presented for Np = 4. This is motivated by on-going studies on the lattice using Kogut-Susskind 
fermions (l6| . 

The running coupling of the Yukawa theory up to O(ay) was computed long ago [34| and reads 

ay (A) - r . (35) 

s+Me log (A /A 

27r 10 § v MS' 



Fig. [S] displays the scale dependence of the Yukawa coupling for different numbers of flavors. This plot illustrates the 
crucial role of Np in the intensity of the interaction as a function of the energy scale and, therefore, in the delimitation 
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of the domain of validity of perturbative calculations. Raising Np from 1 to 6, the energy scale at which the Yukawa 
coupling reaches one increases from ~ 30% to ~ 70% of the maximum scale A^g, a very significant effect. This is 

also clear from Eq. (|35|) . where ay ~ 1/Nf, which strongly suggests that a large Np approach could be appropriate 
for theories with Yukawa- type couplings as has been noticed long ago in a different context by Gross and Neveu |3a| , 
who also identified the need to go beyond the large Np expansion in order to obtain sensible results (see also (l2j|). 
The fermion mass runs according to [l5[ 

m(A) = [ay(A)] 2 < 2JV >+ 3 > m , (36) 

where mo is the value of the mass at the scale where ay = 1, i.e., 

m =m(A ) ; a Y (A ) = l. (37) 

The RG flow of the fermion mass is plotted in Fig. [5] Once again, the number of flavors plays an important 
quantitative role. More interesting, though, is the following feature: in the large Np limit fermion masses tend to 
become scale invariant, i.e., for larger values of Np the running of m becomes flatter. In fact, from Eq. (|36|) . the 
Np — ► oo behavior of m(A) is ~ [(1/Np) log(Ao/A)] 1 / JVF mo, rendering RG corrections to m negligible very fast as 
Np increases. 




The scale dependence of the effective mass of the bosonic field is determined by the following flow equation 

dml N F 2 



SlogfA/A^j 



a Y [m , - 6m ] , 



(38) 



with ay = ay ( A) and m = m(A) given by (|35p and (|36p . respectively. The solution of this equation can be written 
in the following form: 



m|(A) = m 2 (A) 



12iVF 
2A^p + 3 



C 



MS/ 



(39) 



where C is a constant fixed by the boundary condition m^(Ao) = m^.Oi and is given by 

An 



C = 



log 



A 



MS/ J 



2JV F +3 



%,o _ 12N F 
ml 2N F - 3 



(40) 
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It is clear from Eq. (f38|) that the RG flow for exhibits two possible regimes, allowing the boson effective mass 
to increase or to decrease depending on the sign of the term inside the brackets. We illustrate these cases in Fig. 7 
and Fig. 8 for different numbers of fermion flavors. We also show the interpolation between these two regimes in a 
three-dimensional plot, Fig. 9, as we vary the boundary conditions m^fl/rno for Np = 4. For large values of Np, 
the boson mass is strongly modified (cither increasing or decreasing smoothly, depending on the regime) , and the RG 
flow does not seem to become negligible in either case. 




FIG. 10: Pressure versus chemical potential in units 
of the reference scale A^jg- for Nf = 4 for different 
values of A//x in the running of ay. Here, m = and 

?7l0 = 0. 



FIG. 11: Pressure versus chemical potential in units 
of the reference scale A^jg- for Nf = 4 for different 
values of A/fi in the running of ay and m. Here, 
= and m = 0.1Ar-rq. 



Since we consider a general Yukawa theory at finite density, the free parameters are arbitrary and allow for a wide 
spectrum of possibilities for the influence of the RG flow on the equation of state. So, although the RG flow is 
completely well defined by the results presented above, there still remains some freedom in the choice of the following 
scales: A, Ao, Too and m^o- The first usually represents the typical scattering energy scale, and its behavior should 
be motivated physically. In the case of cold and dense matter, the customary natural guess is simply A ~ \x, the 
range of reasonable values for the proportionality constant being determined by physical consistency of the effective 
theory. The others must be fixed by measured observables. Therefore, all of them depend on the specificities of the 
system under consideration. In a given physical system, these parameters can be determined by phenomenological 
constraints, as discussed above, and no ambiguity is left. In the general theory we consider, which describes a whole 
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class of possible physical systems, one can only exemplify the effects of the RG running by choosing a few illustrative 
cases. 

RG flow effects on the pressure are displayed in Figs. 10 and 11 for certain choices of the free parameters, and for 
Np = 4. To simplify the analysis, we opted to keep the bosons massless. In Fig. 10, we show the pressure for different 
values of A//x for massless fermions. It is clear that the pressure is lowered as one increases A/fj,. In the case of 
massive fermions, displayed in Fig. 11, the effect of the RG flow is exactly the opposite. This apparent contradiction 
is not surprising. It is actually a natural consequence of the major influence of the fermion mass on the pressure. The 
fact the the curve with fixed mass m = O.lA-^jg is much lower than the rest is just an artifact of our choice of fixed 
mass being much bigger than the effective mass for these densities. 

V. CONCLUSIONS AND OUTLOOK 

Yukawa theory at vanishing temperature provides at least one of the ingredients for an effective description of the 
thermodynamics of a variety of cold and dense fermionic systems, from condensed matter to particle physics. In this 
paper we have calculated the thermodynamic potential for a general Yukawa theory up to two loops within the MS 
scheme. We have explicitly considered the effect of arbitrary masses for both fields, the spinorial and the scalar, which 
renders the evaluation of two-loop integrals highly non-trivial. Nevertheless, all integrals were computed analytically 
and expressed in terms of well-known special functions. These results were verified against numerical results in several 
different regions of the parameter space. The role of RG running of coupling and masses was also investigated, as 
well as its dependence on the number of fermion flavors. 

Our results show that the pressure is clearly lowered by increasing the mass of the fermions, this effect being 
much more relevant than the ones brought about by variations of the coupling. The effect of the scalar mass is also 
significant, and can be viewed, in a quasi-particle picture, as providing important corrections to the fermion mass 
which, as seen, modifies appreciably the equation of state. The role of the RG flow, on the other hand, is more fuzzy 
in the case of an arbitrary Yukawa theory, and our findings for its influence on the pressure were illustrated in a few 
representative cases. Once the physical system is specified, though, these ambiguities are completely resolved. 

The behavior of running effects with increasing Np is also quite interesting for two reasons. First, it is rather 
convenient for perturbative purposes since ay remains smaller than one even for higher values of the energy scale. 
Then, one could, for instance, compare lattice results to perturbative calculations in a wider range of energy scale. 
Second, because the flow of the fermion mass tends to flatten out, becoming much simpler. Thus, our equation of state 
in the case of Np = 4 can be useful, for instance, in on-going studies on the lattice using Kogut-Susskind fermions 
[lH to deal with effective theories at finite density in the hope of bringing some understanding to the study of the 
Sign Problem [3(| . 

To study the phase structure of systems such as the antiferromagnetic/superconductor in the Hubbard model [l[, 
condensates in the core of neutron stars 0, and QCD, among many others, one has to generalize our description and 
include the presence of a non-zero condensate by computing the full effective potential. This calculation is under way 
for the linear sigma model, and will be presented elsewhere (l8| . 
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APPENDIX: IN-MEDIUM EXCHANGE DIAGRAM WITH MASSIVE FIELDS 

In this appendix we discuss some technical details concerning the explicit calculation of the in-medium exchange 
contribution to the thermodynamic potential with both fermionic and bosonic non-zero masses. 
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1. Matsubara sums 

From Eq. (JSj), solving the trace over Dirac indices and using the following representation of the Kronecker delta 

e/ 3[fc°+p°- M ] _ e /3[p!-f«] 
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Pi - Pi 
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for the bosonic and fermionic Matsubara sums, respectively, with g(q) an analytic function, the triple sum S reads: 
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After a straightforward algebraic manipulation using the properties of the Fcrmi-Dirac and Bosc-Einstcin distribu- 
tions, the above equation can be written as 
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where we used the definitions given in Eqs. (jlip - (fl4|) . 

Finally, the expression in Eq. (JTOJ) is obtained from the result (|A.8[) by taking advantage of the symmetry of the 
integral in Eq. (|A.2|) to conveniently exchange pi <-» P2 in the terms ~ Nf(2). 



2. Renormalization 

Let us now sketch the proof that the UV divergences present in the terms in Eq. (fT0|) that are linear in the statistical 
distributions, 
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belong in fact to vacuum self-energy subdiagrams. 
Using the auxiliary functions [371 ] 
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or, substituting the original expressions (|A.10[) and (|A.ll[t . 
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On the other hand, consider the one-loop diagrams contributing to the amputated bosonic and fcrmionic vacuum 
self-energies. When evaluated on the mass shell they satisfy: 
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Comparing Eqs. (|A.23[) and (|A.24[) with Eqs. (|A.17[) and (|A.18[) , one concludes that L/ and L& are written in 
terms of the vacuum self-energies as stated in Eqs. ([15)1 and (|T5)) . 

Thus the UV renormalization of these terms results directly from adopting renormalized expressions for these 
self-energies. Using dimensional rcgularization in the MS scheme, the measure in momentum integrals is modified as 
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(A.25) 



where e = 4 — rf, with d being the space-time dimension, A the renormalization scale and 7 the Euler constant. A 
standard computation yields for the amputated vacuum one-loop self-energies evaluated on the mass shell the following 
regularized expressions: 
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Defining the integrals over the Feynman parameter x as 
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and substituting the renormalized expressions, obtained from Eqs. (|A.26[) and (|A.27[) by subtracting the poles in 
e = 0, in Eqs. (|15p and (fT6|) . one arrives at the final results in Eqs. (|17|) and (fT8|) . The solution of the integrals on is 
straightforward, yielding Eqs. (fT9|) -([2T |) . 



C. The integral Ji 

Finally, in order to obtain the final analytical result in Eqs. (|2"4]) - (|3l"j) , the integral in Eq. (|23|) must be calculated. 
Explicitly, we have: 
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After evaluating the angular integration, one obtains: 
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where P i = J E 2 — m 2 . 

The integral Xu is simple, yielding: 



In = -\ VPf - m log 
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(33) 
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On the other hand, the integral Tj is extremely more involved. One of the dificulties is the fact that it corresponds 
to a double integral of a log whose argument contains roots which depend on the variable of integration. It is possible, 
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through a convenient transformation of variables, to obtain an integrand without roots. A change of variables with 
this feature was presented in Ref. [24| : 
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Once this new set of variables is adopted, the procedure to solve the integral 2/ includes a sequence of partial 
integrations and several algebraic manipulations until an expression is reached such that it contains only analytically 
solvable integrals, at least for a software of algebraic computation. Solving these integrals, one arrives at: 
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4m 2 



(41) 



where JC\ og is the last integration still to be done 



Kio S (x,z) = ^ \/z(l - z) I dx log 
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(42) 



The solution of /Ci og can be written in terms of products of logarithms whose arguments can assume negative 
values, resulting in a complex result in general. These imaginary parts must cancel out in the end. To verify this 
cancellation and simplify the expression, it is convenient to consider two complementary regimes: z = rn 2 i/Am 2 > 1 
and z = m^/Am 2 < 1, and to use the fact that x = Pf/(n + m) < 1. The final result, after a very long set of 
manipulations, is the one given in Eqs. (1241) (1311) . 
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